Magnetic field induced control of breather dynamics in a single plaquette of 

Josephson junctions. 
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We present a theoretical study of inhomogeneous dynamic (resistive) states in a single plaquette 
consisting of three Josephson junctions. Resonant interactions of such a breather state with elec- 
tromagnetic oscillations manifest themselves by resonant current steps and voltage jumps in the 
current-voltage characteristics. An externally applied magnetic field leads to a variation of the rel- 
ative shift between the Josephson current oscillations of two resistive junctions. By making use of 
the rotation wave approximation analysis and direct numerical simulations we show that this effect 
allows to effectively control the breather instabilities, e. g. to increase (decrease) the height of the 
resonant steps and to suppress the voltage jumps in the current-voltage characteristics. 



Discrete Josephson coupled systems have attracjt 
recently a lot of attention due to the predictior 
and the following observation of intrinsic dynamic lo- 
calized states.u These breather states appear in homo- 
geneously dc current driven Josephson junction ladders 
and arrays in the form of various spatially inhomogeneous 
voltage patterns. These voltage patterns are character- 
ized by a few junctions being in the resistive (whirling) 
state while the rest of all junctions reside in the supercon- 
ducting (libration) states. The breather states manifest 
themselves through additional branches in_jthe current- 
voltage (/-V) characteristics of a systemmD and, more- 
over, can be directly visualized by using the law, tem- 
perature scanning microscopy (LTSM) techniqueJj'B Note 
here, that the origin of such a dynamic localization is not 
the presence of disorder but ajpeculiar interplay between 
nonlinearity and discretenessB, and therefore, breather 
states can play the same important role for underdamped 
Josephson junction arrays as Josephson vortices do for 
parallel arrays a 

The important peculiarity of inhomogeneous resistive 
states is a strong resonant interaction with linear elec- 
tromagnetic excitations, e. g. electromagnetic waves or 
localized electromagnetic oscillations. It was predicted 
in Ref. fj] and carefully studied in Ref. |0 11 that the res- 
onant interaction of the breather state with electromag- 
netic excitations leads to resonant steps (a sharp current 
increase) and various switching phenomena between dif- 
ferent resistive states (voltage jumps) in the I-V curves. 

The simplest case jjthere inhomogeneous resistive 
states can be obtained Jj'EJ is a dc current driven single 
anisotropic plaquette containing three Josephson junc- 
tions, as presented in Fig. 1. It consists of two vertical 
junctions parallel to the bias current 7 direction, and 
a horizontal junction in the transverse direction. The 
dynamics of the system crucially depends on two param- 
eters: the anisotropy 77 = -f^, where I c h and I c y are 
respectively the critical currents of horizontal and verti- 
cal junctions, and the discreteness parameter (normalized 
inductance of the cell), (3l- 
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FIG. 1. Sketch of the plaquette with three Josephson junc- 
tions (marked by crosses) in the presence of uniform dc bias 7 
and an externally applied magnetic field H ex t. Dashed circles 
denote junctions in the resistive (whirling) state. 

In Ref. [ll] we have shown that inhomogeneous resistive 
states can be observed in this system as the anisotropy 
X] is not large. These two breather states are character- 
ized by only one (left or right) vertical junction and the 
horizontal junction being resistive. We found also vari- 
ous further features of the states, e.g. resonant current 
steps and voltage jumps, in the I-V curves. These fea- 
tures can be considered as "fingerprints" of the resonant 
interaction of the breather state with electromagnetic os- 
cillations (EOs) . The voltage positions of the resonances 
depend strongly on the parameters 77 and /3l- 

In this Letter, we report on the properties of breather 
states or even more specifically, on the resonant inter- 
action of the breather states with EOs in the presence 
of an externally applied magnetic field. We show that 
the magnetic field allows to control the strength of the 
resonant interaction just by a fine tuning of the relative 
phase shift between the Josephson current oscillations of 
the two resistive junctions. 

In order to quantitatively characterize the magnetic 
field induced control we write the set of equations for the 
time-dependent Josephson phases < 
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where the nonlinear operator Af(4>) = (f> + cap + sin(0), 
and the dissipation parameter a determines the effective 
damping in the system. Here, we use the inverse of the 
plasma frequency ujj as the unit of time. The externally 
applied magnetic field H ext is characterized by the frus- 
tration / = i- e. the magnetic flux threading the 
cell normalized to the magnetic flux quantum. These dy- 
namic equations have been derived in Ref. [H], using the 
resistivcly shunted model for Josephson junctions, the 
Kirchhoff's current laws and the flux quantization law. 

The Josephson phases are naturally decomposed as fol- 
lows (below we consider a particular breather case as the 
left vertical junction is in the superconducting state): 
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where the breather frequency and the phase shifts c\ 
and c 2 are obtained using a dc analysis^ i. e. neglecting 
the Josephson phase oscillations S(t): 
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In Eq. (||) Sl 2 (t), 8 h (t) correspond to EOs excited in 
the presence of the breather state. By making use of the 
linearization of Eq. (|l|) around the breather state (0) we 
derive two characteristic frequencies of EOs: 
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Thus, as the breather frequency f2 matches any of the 
EOs frequencies, namely Q = w±, a primary resonance 
appears, and the I-V curve displays a resonant current 
step. Using Jie rotation wave approximation method 
elaborated ir£3 we find that the strength of the primary 
resonance is determined by the difference of two Joseph- 
son current oscillations: sin(£lt + c 2 ) — sin(ftt). It leads to 
the magnetic field dependent magnitude of the resonant 
step A7 as 
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Next, we turn to the parametric resonant interaction 
of the breather state with EOs, as £1 = 2cu±. In or- 
der to analyze this case we introduce the new variables 



D{t) = 5 v 2 {t) - S h (t) and S(t) = S$(t) + S h (t). 
equations for D(t) and S(t) are given by: 
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Note-here that we neglected all nonresonant terms in Eq. 

Similarly to a well known case of a parametrically 
driven harmonic oscillatorliJ we find that in the dissi- 
pative case with two degrees of freedom (Eq. (||) ) the 
breather state can become parametrically unstable, when 
the Josephson phases D(t),S(t) ~ e At with A > 0. 
This instability depends on the effective amplitude of the 
parametric Josephson current oscillations h e ff(f) , and 
is determined by the condition: 
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Moreover, the strength of the combination resonance, as 
il = oj + ±w_ , between the breather state and EOs is also 
determined by the parameter h e tf(f). The condition of 
the combination resonance instability is 



h eff(f) > a^/uj+uj- 



(8) 



Thus, the applied magnetic field allows to control the 
strength of the parametric and combination resonances 
through the magnetic field dependent phase shift between 
two Josephson current oscillations, c 2 . 

To test the obtained analysis and to show how the ap- 
plied magnetic field changes the I-V curve we perform 
direct numerical simulations of Eq. (Q). Details of the 
numerical procedure were given elsewhere. EU To estab- 
lish a large dc current bias region where various types 
of resonances occur we use the discreteness parameter 
0l = 1 and the anisotropy 77 = 0.5. We decreased (in- 
creased) the dc bias current and calculated the dc voltage 
V = < 4> 2 >. Numerically simulated I-V curves for 
different values of the applied magnetic field are shown 
in Fig. 2. 

In the absence of the magnetic field ( see Fig. 2a) we 
observed two resonant steps labeled as A and B, and the 
switching phenomenon (the voltage jump to the homoge- 
neous resistive state) labeled as C. All types of resonant 
interactions of the breather state with EOs are present 
in this I-V curve: the parametric resonance A leads to 
the the resonant step at upper values of the dc bias 7, 
the primary resonance B leads to the resonant step at 
lower values of the dc bias 7, and the voltage jrump C is 
the consequence of the combination resonance. EJ 
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FIG. 2. I-V characteristics of a breather state for /3l = 1, 
a — 0.1 and 77 = 0.5. The values of the magnetic field are (a) 
/ = 0; (b) / = -0.1; (c) / = 0.4. The dotted lines show the 
dependence of the characteristic frequency combinations of 
EOs on the dc bias 7 (Eq. (0) ). Arrows indicate the various 
switching processes. 

In the presence of an applied magnetic field we find 
that the I-V curve changes drastically. Thus, for a small 
magnetic field / = —0.1 the lower voltage resonant step 
B practically disappears, but the switching phenomenon 
C becomes even stronger (Fig. 2b) . In the opposite limit 
of a large magnetic field, as / = 0.4, the parametric and 



the combination resonances (A, C) disappear from the I- 
V curve but the magnitude of the lower voltage resonant 
step B (primary resonance) increases (Fig. 2c). 

The dependence of the magnitude of the lower resonant 
step on the frustration displays a minimum on small neg- 
ative values of / and reaches a maximum as / ~ 0.4 
(circles, in Fig. 3a). We find that the outcome of our nu- 
merical simulation is in good accord with the theoretical 
analysis (Eq. (JsJ) and solid line in Fig. 3a). 




FIG. 3. The dependence of the magnitude of the resonant 
steps on the frustration /: (a) the lower resonant step (pri- 
mary resonance, B). The numerical results (circles) and the 
analytical prediction (Eq. (^) ) are shown correspondingly by 
circles and solid line, (b) the upper resonant step (parametric 
resonance, A), the numerical results (circles) and the analyti- 
cal prediction (solid line) based on the expression for h e ff(f) 
(Eq. (Q) ) are shown. 

The numerically simulated dependence of the mag- 
nitude of the upper resonant step A is shown in Fig. 
3b. This dependence displays a peculiar oscillation and 
is in qualitative agreement with the magnetic field de- 
pendence of h e ff(f). Moreover, as suggested by the in- 
equality in Eq. (pj) is valid, in the region of frustration 
0.2 </ < 0.37 a specific "window" of magnetic field oc- 
curs. In this region of the magnetic field the parametric 
resonant interaction is very weak and the resonant step A 
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disappears from the I-V curve. Note here that the width 
of the "windows" increases with the damping parameter 
a (Eq. (0) ) , and this allows for an accurate measurement 
of damping at low temperatures as a is very small. 

In the same region of frustration, a similar "window" 
exists for the combination resonance and correspondingly 
the switching of the breather state to the homogeneous 
resistive state is suppressed. This effect has a peculiar 
consequence. Because of the instability of the breather 
state and the corresponding switching to the homoge- 
neous resistive state ( see Fig. 2b), it can be very hard 
to experimentally observe the lower part of the I-V curve. 
However, in the presence of a magnetic field the breather 
can pass this region of dc bias current and penetrate to 
the "forbidden" current region. 

In conclusion we presented an analysis and direct nu- 
merical simulations of the inhomogeneous resistive state 
in a single plaquette consisting of three Josephson junc- 
tions. We find that the I-V curve of such a breather 
state displays peculiar features, namely resonant steps 
and voltage jumps. The magnitudes of all features are 
determined by the phase shift between two Josephson 
current oscillations that in turn, can be controlled by an 
externally applied magnetic field. 
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